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Abstract 

A system of two parallel Josephson junction arrays coupled by interlayer 
capacitances is considered in the situation where one layer is in the vortex- 
dominated and the other in the charge-dominated regime. This system shows 
a symmetry (duality) of the relevant degrees of freedom, i.e. the vortices 
in one layer and the charges in the other. In contrast to single-layer arrays 
both contribute to the kinetic energy. The charges feel the magnetic field 
created by vortices, and, vice versa, the vortices feel a gauge field created by 
charges. For long-range interaction of the charges the system exhibits two 
Berezinskii-Kosterlitz-Thouless transitions, one for vortices and another one 
for charges. The interlayer capacitance suppresses both transition temper- 
atures. The charge-unbinding transition is suppressed already for relatively 
weak coupling, while the vortex- unbinding transition is more robust. The shift 
of the transition temperature for vortices is calculated in the quasi-classical 
approximation for arbitrary relations between the capacitances (both weak 
and strong coupling). 

Two-dimensional (2D) Josephson junction arrays have attracted much attention because 
of the experimental and technological progress and the rich underlying physics (see Ref. |I| 
for review). Classical 2D Josephson junction arrays, where the Josephson coupling energy Ej 
between the superconducting islands dominates, is a standard example of a system exhibiting 
the Berezinskii-Kosterlitz-Thouless (BKT) transition - the unbinding of vortex-antivortex 
pairs at a certain temperature Tj PJ^. The transition separates a superconducting phase 
at T < Tj (X Ej, where vortices are bound, from a resistive phase. It was realized later (see 
e.g. Refs. 0-0]) that charging effects, associated with the capacitances of the islands to the 
ground Cq and of the junctions C, lead to quantum fluctuations of the phase and suppress 



the BKT transition temperature. Beyond a critical value of the charging energy Ec = 
min{e^/2C, e^/2Co} the transition temperature vanishes, and the superconducting phase 
ceases to exist. The next step was to understand that in the "extreme" quantum limit 

Ec ^ Ej, where the quantum fluctuations of the phases are very strong, the vortices are 
ill-defined objects. In this regime the charges on the islands become the relevant variables. 
If, furthermore, Co ^ C, the interaction between the charges is (nearly) the same as that 
of the vortices in the quasi-classical array. In particular, the charges can be considered as a 
2D Coulomb gas |T^, and they undergo a BKT transition at temperature Tq oc Ec- The 
phase below the transition is insulating. A finite value of the Josephson coupling between 
the islands suppresses this transition. As another example we mention the influence of 
dissipation (e.g. Ohmic dissipation) on the phase transition in the array, which was first 
noted in Ref. |T2|. We are not going to review these theoretical results, however it is necessary 
to stress that the theory of 2D Josephson junction arrays is far from being settled. 

Below we describe another, more complicated system - two parallel 2D Josephson junc- 
tion arrays with capacitive coupling between them (no Josephson coupling^). Probably the 
most interesting situation arises when one array is in the quasi-classical (vortex) regime 
while another one is in the quantum (charge) regime. Then the vortices in one layer and the 
charges in the other one are well-defined dynamical variables. Another important feature 
of the present system is that the strength of interaction between charges and vortices is 
controlled by the interlayer coupling and consequently this interaction may be weak or 
strong, whereas in usual 2D Josephson junction arrays the strength of charge- vortex inter- 
action is of the same order as either the charge-charge or the vortex-vortex interaction. We 
also show that the physical realization of this interaction is rather different from that in one 
array. Hence, at least for weak interlayer coupling, one should expect two BKT transitions, 
the first for charges in one layer, and the second for vortices in the other one. In this article 
we provide the theoretical description of the coupled system and calculate the shifts of the 
transition temperatures due to the interlayer interaction. 

We consider two parallel Josephson junction arrays, i.e. (square) lattices of supercon- 
ducting islands connected by Josephson links. As usual, we suppose that the magnitude 
of the order parameter in the islands is constant while its phase fluctuates from island to 
island. The partition function of the system may be expressed conveniently in terms of these 
phases (the indices i label the islands in each array and /i = 1, 2 refers to the number of 
array) 

Z = \{p d<f>^^d<f>§ J DMr)DMr) exp(-5{0}) . (1) 

« ^ {mii},{mi2} 

Here the path integration over phases is carried out with the boundary conditions 

MO) = <pfj; MP) = <Piu + 2vrmi^, 



Multi-layered systems with Josephson coupling between layers have been discussed in the liter- 
ature (see e.g. ||T^ , [l5| ). The analog of the BKT transition is in this case the disruption of vortex 
rings. In the limit of weak Josephson couplings this system is reduced to the 2D XY-model, while 
in the opposite case of strong coupling it is essentially the 3D XY-model. This situation, however, 
is absolutely different from the one we describe below. 
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where (3 is the inverse temperature, and h = 1. The Euchdean effective action S{(j)} is 



+ ^ E('^il - ^^^y + Ejl E[l - COS(0.1 - 0,1)] + ^J2 Efl - COS(0.2 - 0,2)] ^ . (2) 



Here Cq^ are the capacitances of the islands in the array /i relative to the ground, are 
the capacitances of the junctions in the array fi, and Cx are the interlayer capacitances. 
Furthermore, Ej^ are the Josephson coupling constants in the layers. Here and below we 
use the symbol to denote the summation over nearest neighbors only, and each pair is 
counted once; the symbol stands for the summation for all values i and j (in particular, 
each pair except {ii) is counted twice). 

From now on we choose parameters such that the array 1 is in the charge (quantum) 
regime while the array 2 is in the quasi-classical (vortex) regime. In terms of the phase 
variables this means that in the array 1 the phases on each grain are strongly fluctuating 
in time, while in array 2 they are nearly time-independent. This regime is described by the 
conditions 

Eji « eye,, Ej2 » eyC2, 

with = ma:x{Co^,C^,Cri.}. Below we first calculate the shift of the BKT transition 
temperature for vortices in the array 2. This does not require the introduction of charges 
and vortices and may be done in the phase representation. Then, we turn to the BKT 
transition for charges in the array 1. For this purpose we move from a description in terms 
of phases to one in terms of charges and vortices, and use the duality of the resulting action 
to investigate the transition. At the same time, we will show that charges and vortices in 
this system can be considered as two-dimensional dynamical particles with masses. The 
charge-charge and vortex-vortex interaction are essentially those of 2D Coulomb particles, 
while the charge-vortex interaction is more peculiar. 



BKT transition for vortices 

The shift of the BKT transition temperature for vortices in the array 2 due to the 
coupling to the array 1 can be calculated easily if we set the small parameter Eji to zeroQ. 
Then the action for the phases (pn becomes Gaussian and the latter may be integrated out. 
After that the shift of the BKT temperature for vortices may be obtained by means of the 
quasi-classical expansion []7|,p!6 . 



The first step requires a comment. The path integration over the phases of the islands 
()ji(r) in the array 1 is, as usual, performed by a linear shift of variables in order to eliminate 



^This means, in particular, that the results obtained below are valid also in the case when the 
array 1 is in the normal state. Because of the e-periodicity in this case the boundary conditions in 
the array 1 should read 4'ii{(3) — 0ji(O) = Airmn. As we will show, this does not change the final 
result. 
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terms linear in However, the new (shifted) variables do not generally satisfy the bound- 
ary conditions, and consequently the integration is not possible. If, nevertheless, the array 2 
is in the quasi-classical regime, the contributions of all non-zero winding numbers to the 
partition function are exponentially small in comparison with the contribution of = 
(see, e.g. If we neglect these small contributions, the phases 0i2 become periodic, and 

the boundary conditions in the array 1 are met automatically. After integration over the 
phases (pn we find 

^ = n fj J DMr) exp(-S{0,2}) (3) 
with the effective action 

S{4>} = l^dTl-^Y. MMjj2 + Ej, ^[1 - COS(0,2 - <Pj2)] \ . (4) 
° I ^ U (ij) J 

Here [A2]ij is the effective capacitance matrix for the layer 2 (see also below) 

[^2h = -^^[Q^% + C,[Q,],, (5) 
and the matrices (/i = 1,2) have a form 



[Q 



4 + i = j 

— 1 i and j are nearest neighbors 

otherwise 



Since the array 2 is supposed to be in the quasi-classical regime, only weakly time- 
dependent periodic paths 4'i2{'^) are important. Hence we may write the phases in the form 
0i2(T) = + /i(r), where 

oo 

/i(r) = (3~^ "^[fi^Un) exp{-iuJnT) + f*{uJn)exp{iUnT)]. 

n=l 

is expressed as a sum over Matsubara frequencies ci;„ = 27mj3^^. Now the action may be 
expanded to quadratic terms in /i(r), yielding 

^{^S?, /.(r)} = PEj, ^[1 - cos(4°^ - 4?)] + 

(ij) 



+ ['drl^Y: m^l^J, + ^ UMr) - fAr)r cos(4°) - ^J?) I . (6) 
JO I »e ■■ z ^..^ I 



Note that the first term is the classical action of 2D Coulomb gas . Finally, one performs 
the cumulant expansion |?,16] in the last term in brackets in Eq. (^. As a result the action 
has exactly the same form as the classical one, but with the renormalized temperature 

(3-^r'-U'dr{{Mr)-f,{r)f). (7) 
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Here the angular brackets denote the averaging with the effective action 

ij 

From this we obtain the shift of the transition temperature due to the interlayer coupling 

Tj = Tjo - e^A/3, (8) 

where 

and the second term is the matrix element taken for neighboring islands. (We have assumed 
that the matrix A2 depends on the distance between the islands only). Here Tjq is the tran- 
sition temperature for a classical 2D Josephson junction array |^ (to be of order Ej2). Eq. 
(§) is the result we were aiming at, however in order to obtain some analytical expressions 
we evaluate the quantity A in some approximations. In Fourier representation the matrices 
Qfj,{k) have the form 

g^(k) = (kay + iC,)-\C, + Co^), ka < 1, 

with a being the lattice parameter. Consequently the matrix A^^ is 

\-Uk] = Cijka)^ + + Cqi 

' ^ ' (Ci(A;a)2 + a + Coi){C2{kay + + C02) - ^ ^ 

If we replace the first Brillouin zone by a circle cut-off at A; < a~^, the integration over the 
angular variable can be performed easily, and we obtain 

A = — (^'\dkk2^{k)\l - Jo(H]- 
27r JO 

In the range of integration the Bessel function Jo can be approximated by its expansion 

Jo(x) ^ 1 - 

Finally, in the case Cq^ ^ (this situation is the most interesting) we obtain 

j,^ _rp^ _ 6^ 1 _ — In C'xC'i + CxC2 ^^Q^ 

487rCe//' Ceff C2 C2 CxCi + CxC2 + C1C2 



It is seen that the effect of layer 1 is merely the renormalization of the effective capaci- 
tance. As a result, the BKT transition for vortices in the layer 2 is suppressed. We should 
emphasize that the result (|l3) is valid for arbitrary capacitances C^, Ci and C2. The only 
restriction is the validity of the quasi-classical approximation. The shift of the transition 
temperature should be small, or, in other words, e^/Ce// <^ Ej2. In particular, for ^ C2 
(weak coupling between the layers) one obtains Ceff ~ C2 irrespectively of Ci — vortices 
in layer 2 do not feel the presence of the layer 1. In the case Ci = C2 -C Cx the effective 
capacitance is C^ff = 2C2, while for Ci = ^ C2 one has Ceff = 2Cx/3. It is seen that 
in the latter case the temperature begins to feel the presence of the first layer, however the 
absolute value of the shift becomes now small. 
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Charge — vortex duality and BKT transition for charges 



Before we turn to the description of the BKT transition for charges, it is necessary to 
stress the following. As shown by the effective action (|^), the interlayer capacitance Cx not 
only couples the layers, but also renormalizes the capacitances Coi and C02 of the islands 
to the ground. Hence the logarithmic interaction between the charges in each layer has a 
finite range for any non-zero due to the screening, and the BKT transition, is, strictly 
speaking, absent. One should realize, however, that in the situation Cqi ^ Ci the 

screening length ^1 ~ a{Ci/Cxy^^ can be very large. Below we assume that these inequalities 
are satisfied and the range of interaction ,^1 is large enough to make it meaningful to speak 
about the charge-unbinding transition. (Note that like any phase transition in a finite system 
this transition is smeared; in other words, the resistance grows exponentially, and, strictly 
speaking, for finite ^ it is impossible to distinguish between insulating and resistive phases). 
On the other hand, already for relatively weak coupling ~ Ccr ^ Ci this description 
becomes meaningless and the insulating phase is absent. As we show below, in the small 
range Cx ^ Ccr the transition temperature for charges in layer 1 does not feel the presence 
of the layer 2 and hence is essentially the charge-unbinding temperature for one Josephson 



junction array [jTO[. Nevertheless, the charge- vortex description required to obtain this result 
gives rise to an interesting physical model to be described below. 

Now we move from the phase description (|1]),(0) to a charge- vortex description. First we 
introduce the large capacitance matrix 

Here is the capacitance matrix in the array fi while Cx = CxSij. The inverse matrix in 
the Fourier representation reads as 

f-'^Clr] = i X ('I2'l 

(C^{kay + Cx + Coi)iC2ikay + Cx + Co2)-C^, ^ ^ 



C2{kay + Cx + C02 Cx \ / Ar^(A;) A-\k) 

Cx C^{kaf + Cx + C,^)-\ k-\k) A2\k) 



Indices /x, z/ = 1, 2 again label the array. The matrix C~^ describes the interaction of charges. 
We have also introduced for later convenience the matrices A~^, describing the interaction 
of charges within layer n, as well as A~^ referring to charges in different layers (cf. Eq. (^. 
Then the effective action (0) can be rewritten in terms of integer charges g^^ of each island 
and phases (pi^ 



-^0 I ■;,,,/ 



+ 



+Eji cos(0ii - 0,1)] + Ej2 ^[1- cos(0i2 - 0,2)] [ . (13) 

Now it is possible to introduce vortex degrees of freedom by means of the Villain transfor- 
mation [|l^] (see also |jl8|). It is important that this procedure deals only with the phase 



6 



variables and does not affect the charge interaction (the first term in Eq. ([T3D). The phase 
terms (the second and the third one in Eq. ([T3|) ), however, have exactly the same form 
in the problem of two arrays as for a single Josephson junction layer. The procedure for 
a single-layer array array is discussed in details in Refs. [p!0| , |TT| , the generalization to the 



double-layer system is straightforward. The partition function becomes 

^ = n E E exp(-5{g,4), (14) 

* {qii,qi2} {i'ii,Vi2} 

where the effective action for integer charges g^^ and vorticities Vi^ is 

SU, v} = [ dr Ue'Y.T. ?^.(^) \C-X^A^) + . F F ^ E Ur)Gf;>U^) + 

+ r4^ E <i^2ir)G\fq,2{r) + nEj,F{e,Ej,) ^ VaG^^^v,, + nEj^ ^ v,2G\fv,2+ (15) 

+T' E + E ^i2{'r)^ijVj2{ 

ij ij 

Here we introduced the discrete time variable; the time lattice spacing in the array is 
of order ~ (SEj^Ec^)^^^"^ , Ec^ = e^/2C^. The time integration and derivatives are 
continuous notations for a summation over time lattice and for a discrete derivative 

respectively. The function 

^^^^ " 2xln(Jo(x)/Ji(x)) ^ 2xln(4/x)' ^' 

is introduced to "correct" the Villain transformation for small Ej [|r^ . As we see, its entire 
effect is to renormalize (to increase) the Josephson coupling in the layer 1; the renormalized 
coupling Eji reads as 

Ej, ~ [SEnEciY^' {ln{Ec,/En))-' . (16) 

Note that for Eji <^ Eci (charge regime) one obtains Eji ^ Eci- 
The kernel 

Vi - Vj N 



0jj = arctan( 



Ob 2 Ob -i 



describes the phase configuration at site i around a vortex at site j. Finally, the kernel G^^^ 
is the lattice Green's function, i.e. the Fourier transform of At large distances between 
the sites i and j it depends only on the distance r between the sites and has a form (see e.g. 

G;f = ln(e», a«r«e^, i, = a{CJG^^/\ (17) 

Later on, we assume that the linear size of each array is much less that the range of interaction 
This means, in particular, that we assume G^ -C 6*2- 
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The action (|T5D depends on the charges and vorticities in both layers. However, in our 
situation, when the layers 1 and 2 are in the charge and vortex regimes, respectively, the 
vortices in the layer 1 and the charges in the layer 2 may be integrated out . To do this we 
suppose the latter variables to be continuous (strongly fluctuating), and neglect the kinetic 
term for charges in the layer 2 (g2C^^-'g2)- Then after performing the Gaussian integration 
we obtain the effective action for charges qn in the layer 1 and vorticities Vi2 in the layer 2 
(to be referred below as qi and Vi) 



TT 



8E, 



C2 



E 



ViiT 



1-j 



Att^GiG2 



kl 



1 ijk 



To derive Eq. (p!8[) we have taken into account the explicit expression for the large capaci- 
tance matrix (P^. 

The action (p!8D is the central result of this section. It looks rather similar to the effective 
charge- vortex action in one Josephson junction, but the most important difference is that 
while in one layer either charges or vortices are well-defined, 
of well- defined dynamic variables on each site 



Eq. 



' p!8D describes the system 
charges in the layer 1 and vortices in the 
layer 2. We postpone the discussion of physics in this system until the next section, however 
it is clear that the action shows a duality between charges and vortices. The second term 
in the square brackets is small if Cx -C Ci, C2. Both kinetic terms for charges and vortices 
violate the duality due to the numerical coefficients. However close enough to the transitions 
these terms produce only small renormalization of the transition temperature, and are not 
important. Another interesting feature of this action is that the last term, describing the 
interaction between charges and vortices, is also small, while in a single-layer array the 
interaction is always of the same order of magnitude as another terms. 

It is obvious that for long-range interaction of the charges in the layer 1 they also ex- 
hibit the BKT transition, and under the conditions where the action (p!8|) was obtained the 
transition temperature does not feel the presence of the layer 2: 



CO 



c 



Eji 
24^' 



Charge and vortex motion 

To understand the physics described by the action ( [T8| ) it is instructive to map this 
model onto the 2D Coulomb gas. For this purpose we move from the space-time lattice to 
the continuous medium and introduce the coordinates of the vortex centers and charges 

m 

v^{r)^Y^nS{r-R4T)). (19) 

n 
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Here qm = ±1 and f„ = ±1 represent charges and vortices respectively; Tmir) and Rn(r) are 
the corresponding coordinates of the charges and of the vortex centers. Now the partition 
function reads 



oo oo „ 

^= E E Dv,{T)...Dv2M{r)DRi{T)...DR2N{r)exp{-S{r,R}), (20) 

A/=0 Af=0 

and we are going to deal with the effective action ^jr, R}, describing the behavior of the 
system of 2M charges (of which M are positive, q = 1, and the other M are negative, 
q = —1), and of N positive {v = 1) and N negative {v = —1) vortices. 

The first and third terms of the action (ll8|) can be easily transformed by means of 
decomposition (O). The first one produces the potential energy of charge interaction. 



S^it = Tdr f: qmqnG^'HrUr) - r„(r)). (21) 

^ ^0 m,n=l 

In principle, the summation includes the terms with m = n; these, however, may be excluded 
from this sum, giving rise to the chemical potential for charges. The third term in Eq. (|18D 
yields the interaction of vortices 

Stl = 7rEjJ dr v^v^G'''\RUr)-Kn{r)). (22) 
Jo 



m,n=l 



Here again the term with m = n gives rise to the chemical potential for vortices. The terms 
(^) and (p2|) are essentially the action for (classical) Coulomb gases of charges and vortices, 
respectively 0. 

If we neglect the small correction proportional to the CI/C1C2 in the fourth term in 
Eq.(P!HD then the second and fourth terms can be transformed to the kinetic energy of 
charges and vortices respectively |jlO|- The second term gives 



1 rP 2M 

(23) 

Ejl "'0 m,n=l 



2M 

oil) 



We have introduced the mass tensor 19 



M^5(r) = -V^VsG^^\r). (24) 

It decreases proportional to for r ^ a, and consequently may be approximated by a 
local function 

Then the kinetic term for charges takes a simple form 

2M 



kin - ^ / E ^m(^)- (25) 



c{q) 

2a^Eji Jo 



m=l 



Similarly, the fourth term in Eq.(0) produces the kinetic term for vortices 
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C2 



I drY^Rlir). (26) 



Finally, the last term in Eq.(|T8D) is responsible for the interaction between charges and 
vortices. The corresponding term in S{r,R} is 

S,. = -J^ [drY^v^qn I rfr'VR^e(R„ - r')G(^)(r' - r„)R™(r). (27) 

The integral over r can be calculated explicitly, yielding 

Sqv = - rfr^ w„R^(r)A„(R„) (28) 
Jo „ 

with 

Am(Rm) = J2 ^n^i^rnir) - r„(r)), 

n 

a(r) = -^^(l + 21n(C,/r))[zxr]. 

The resulting action is 

S{r, R} = S^^f. + + Slf^ + Sl^l^ + S'^t,. (29) 

The action ( pQ]) is essentially that of two 2D Coulomb systems. The charges and the 
vortices can be viewed as particles with masses 

respectively. Charges interact via the effective capacitance, vortices via the usual logarithmic 
interaction with strength Ej2. Furthermore, the vortices produce the vector potential a for 
the charge^; the magnetic field associated with this vector potential is 

Its sign depends of the signs of the corresponding vortex and charge. Apart from its quite 
peculiar functional form, another important feature of this field is the small factor C^/Ci. 



^This seeming asymmetry is rather artificial. In Eq.([l8|) one can rewrite after a partial integration 
the charge- vortex interaction term in order to obtain the vector potential for vortices, created by 
charges, as well. However, this vector potential contains always the small factor Cx/Ci. 



10 



Summary 



We have investigated the system of two 2D Josephson junction arrays coupled by capac- 
itances Cxi in the situation when the arrays 1 and 2 are in the charge and vortex regime, 
respectively. In the case of weak coupling C-^ <S Ci, C2 the system shows an (approximate) 
duality between dynamical charges in one layer and dynamical vortices in the other one. In 
contrast to a single layer array, both variables are well-defined. The system is equivalent to 
two 2D Coulomb gases of massful particles. The charges feel the magnetic field created by 
vortices, and, vice versa, the vortices feel the gauge field created by charges. In this respect 
the system resembles the composite fermion model of the fractional quantum Hall effect, 
however the magnetic field is now small and has another functional form, so one may expect 
different physics. In this regime the system shows two BKT transitions, one for charges 
and another for vortices, and the coupling between the layers suppresses both transitions. 
Although one could expect the suppression of one transition and the enhancement of an- 
other one, the suppression of both transitions is rather natural, since the capacitance Cx also 
renormalizes the capacitances of the islands to the ground. The BKT transition for charges 
vanishes even for very small values of Cx-, however, the BKT transition for vortices survives 
under condition /Ceff <C Ej2 irrespective of the relations between the capacitances Ci, 
C2 and Cx- The shift of this temperature due to the capacitance effects is calculated within 
the phase representation for both cases of weak and strong coupling. The effect of the layer 
1 is to renormalize the capacitance matrix in the layer 2. For weak coupling C^ <^ C2 (irre- 
spective of Ci) the vortices do not feel the presence of another layer, and the temperature 
remains the same as for one layer. However, for Cx S> C2 different situations are possible. 

In summary, we would like to emphasize that the system of two coupled Josephson 
junction arrays may exhibit quite rich and interesting physics. We have investigated some 
limiting cases, however, the further rich behavior of this system can be expected in other 
cases. In particular, the magnetic field created by vortices seems to be rather unusual and 
interesting. We hope that experimental studies of this system will be performed in the near 
future. 
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